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Abstract. A second neighbor t— J+V model for the chain subsystem of the Sr\4,Cu2i04,\ has been extracted 
from ab-initio calculations. This model does not use periodic approximation but describes the entire chain 
through the use of the four-dimensional crystallographic description. Second neighbors interactions are 
found to be of same order than the first neighbors ones. The computed values of the second neighbors 
magnetic interaction are coherent with experimental estimations of the intra-dimer magnetic interactions, 
even if slightly smaller. The reasons of this underestimation are detailed. The computed model allowed 
us to understand the origin of the chain dimerisation and predicts correctly the relative occurrence of 
dimers and free spins. The orbitals respectively supporting the magnetic electrons and the holes have been 
found to be essentially supported by the copper 3d orbitals (spins) and the surrounding oxygen 2p orbitals 
(holes), thus giving a strong footing to the existence of Zhang- Rice singlets. 



PACS. 7 l.lO.Fd, 71.27.-ha, 71. 23. Ft 



1 Introduction 



One-dimensional quantum systems have attracted a lot 
, of attention in the past decade due to the large diversity 
' of their low energy physics. In particular, spin-chains and 
spin-ladders systems have been extensively studied. The 
characteristic of the Sri±- x A x Ou,2iOn {A = Ca, Ba, Y, Bi, 
etc.) family of transition-metal oxides is that they are com- 
posed of both spin-chains and spin-ladders sub-systems. 
These compounds are formed of alternated layers (in the 
(a, c) plane) of each of the two subsystems p^. Both lad- 
ders and chains are in the c direction, however their re- 
spective translation vector (c c and ci) are incommensu- 
rate. In the pure compound, Sr 14CU24O41, the layers are 
largely separated (~ 3.3 A) and considered as electroni- 
cally non interacting. Nevertheless, the low energy proper- 
ties of St\/JJu2aP^\ agree neither with those of spin chains, 
nor with the properties of spin ladders. 

St\^Cu2aO/±\ is a semiconductor with a 0.18 eV gap |2] 
at T < T* — 250 K. The spin ladders have a singlet 
ground state with a spin gap of about 35 — 47meV j3] 
Surprisingly the spin chains also exhibit a singlet 
ground state with a spin gap of 11 — 12meV P]|^]|rj][7||^j 
while homogeneous spin chains are known to be gap-less in 
the spin channel. Susceptibility and ESR measurement 
E] suggested that the spin gap in the chains is due to the 
formation of weakly interacting spin dimers. Neutron scat- 
tering experiments [3] have later confirmed their existence. 



The spins are supported in the chain subsystem by the 
3d ac orbitals of the Cu 2+ ions, and in the ladder subsystem 
by the 3d a 2_ c 2 orbitals of the Cu 2+ ions. In the chain sub- 
system the magnetic orbitals are coupled via two nearly 
90° Cu-O-Cu bonds. Let us note at this point that, in such 
geometries, the super-exchange paths through the oxy- 
gen orbitals interact destructively and therefore nearest- 
neighbor (NN) exchange interactions are expected to be 
small and ferromagnetic. In the ladder subsystem, the pic- 
ture of the NN interactions is very different since there are 
mediated via nearly 180° Cu-O-Cu angles. Such geome- 
tries are known to produce strong super-exchange mech- 
anism via the bridging ligands and thus large anti-ferro- 
magnetic interactions. 

Formal charge analysis shows that the Sri4C^2404i 
compound is intrinsically doped with six holes by for- 
mula unit (f.u.). Similar to high-T c superconductors, the 
holes were expected to be mainly supported by the oxygen 
2p orbitals and to form Zhang-Rice ^U] singlets with the 
associated-copper hole. NEXAFS (Near Edge X-ray Ab- 
sorption Fine Structure) experiments |11) later supported 
this assumption. 

It has been established, from neutron scattering [3] 
and X ray spectroscopy ^2] experiments, that the chain 
dimcric units are formed by two next-nearest-neighbor 
(NNN) spins separated by a hole. Cu NMR measurements 
exhibited the presence of two kinds of holes on the chains 4 , 
namely with intra- and inter-dimer localization. The rela- 
tive occurrence of the two types of holes (0.65 ~ 1/2) led 
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the authors to propose a charge-order model with dimers 
separated by two holes. This assumption has been con- 
firmed by neutron scattering experiments |7lll3j that have 
shown to be consistent with a five units periodicity. Such 
a picture leads to a chain filling of 6 holes per f.u., that is 
with all the holes located on the chains. 

The question of the holes repartition between the chain 
and ladder subsystems is however still under debate. In- 
deed, NEXAFS experiments ^I] evaluated to 0.8 the num- 
ber of holes on the ladders legs oxygens. Magnetic suscep- 
tibility measurements %■ exhibit a filling of 3.5 spins per 
f.u., that is 0.5 more holes than the maximum number 
given by the formal charge analysis. Finally, magnon-hole 
scattering experiments ^3] suggest a total localization of 
the holes on the chain subsystem for temperatures lower 
than T* ~ 200 K. 

The origin of the chain electronic dimerisation has 
only recently being elucidated. The hypothesis of a spin- 
Peierls transition has been rapidly eliminated since the ex- 
pected signatures in the magnetic susceptibility and spe- 
cific heat were not found. It thus has long been supposed 
that the dimerisation originates from the competition be- 
tween first and second-neighbor spin interactions as pre- 
dicted by Majumdar and Gosh [23 and Haldane [T^]. We 
recently showed, using ab initio calculation, that the origin 
of the dimerized state is of structural origin, even if not of 
spin-Peierls one J7| . Indeed, this is the structural incom- 
mensurate modulations of the chain subsystem, with the 
periodicity of the ladder one, that strongly modulate the 
spin/hole orbital energy and localize the spins in a dimer 
pattern. 

These incommensurate structural modulations, that 
are most of the time neglected, have thus proved to be 
crucial for the low energy properties of this system. Inde- 
pendently to the on-site orbital energy modulations, the 
structural distortions of the chain subsystem were shown 
to be responsible for strong modulation of the first neigh- 
bor effective interactions \Tj\, namely NN hopping and 
NN exchange integrals. One can thus suppose that these 
structural modulations will also affect the second neighbor 
interactions that are expected to be of the same order of 
magnitude than the NN ones in these 90° oxygen bridged 
systems. Indeed, while the super-exchange path through 
the bridging oxygens is hindered in the NN interactions, 
this is not the case for NNN ones where interaction paths 
through two oxygen bridges are favored. In addition, these 
second neighbor effective integrals are crucial for the low 
energy properties of this system since the structural mod- 
ulations localize the spins into NNN dimers. 

The present work thus aims at giving a complete de- 
scription of the chain subsystem within a second neigh- 
bor t — J + V model, taking explicitly into account the 
structural incommensurate modulations. The incommen- 
surate character of the problem will be treated within 
the four-dimensional representation and we will describe 
the electronic structure without the use of a periodic ap- 
proximation. Using this complete model, we will further 
analyze our results as a function of the four crystallo- 
graphic coordinate and in particular study the different 



(quasi)independent spin entities, their relative proportion 
and positions. This four dimensional method presents the 
advantage of giving a complete description of these in- 
commensurate systems, without the use of any periodic 
approximation. 

The next section will describe the ab-initio methods 
used in the present calculations. Section III will analyze 
the nature of the holes and yield ab-initio evidence of 
Zhang-Rice singlets. Section IV will give the second neigh- 
bor t — J + V model as a function of the fourth crystal- 
lographic coordinate r. Section V will details the NNN 
double-bridged super-exchange mechanism and its calcu- 
lations. Section VI will analyze the exchange fluctuations 
along the chains. Finally the last section will conclude. 

2 The ab-initio method 

It is well known that magnetic and transfer interactions 
are essentially local in strongly correlated systems and 
can thus be accurately evaluated using embedded frag- 
ment ab-initio spectroscopy methods ^H]. The long range 
electrostatic effects are treated within a bath composed 
of total-ion pseudo-potentials an d charges. The open- 
shell character of the magnetic/hole orbitals, the strong 
electronic correlation as well as the screening effects are 
efficiently treated using quantum-chemistry ab-initio spec- 
troscopy methods. The present calculations were performed 
using the Difference-Dedicated Configuration Interaction 
method 20 (DDCI) . This configuration interaction method 
has been specifically developed in order to accurately cal- 
culate the low energy states (spectrum and wave func- 
tions) of open-shell, strongly correlated formally finite sys- 
tems. Used on embedded crystal fragments, the DDCI 
method allow to obtain effective local interactions with a 
great accuracy and reliability. For instance, it was proved 
to be very efficient on copper and vanadium compounds 
such as high T c copper oxides [23 or the famous a'NaSfaOs 
compound j22| . 

The quantum fragments are defined so that to include 
(i) the magnetic centers, (ii) the bridging oxygens me- 
diating the interactions, and (iii) the first coordination 
shell of the preceding atoms which is responsible for the 
essential part of the screening effects. Second-neighbor in- 
teractions are thus extracted using three centers Cu^Og, 
fragments (see figure QJa). NNN exchange interactions J2 
are obtained from the doublet-quartet excitations energies 
when three magnetic electrons are considered. NNN hop- 
ping and first neighbor bi-electronic Coulomb repulsion 
are obtained from the 3 singlets and 3 triplets of the frag- 
ment containing one magnetic electron less. Let us point 
out that the 3 centers calculations also yield the NN inter- 
actions, which were previously obtained using two centers 
C112O6 fragments |17|. The comparison between the eval- 
uations of the first-neighbor integrals obtained from the 
2 centers and 3 centers fragments allows us to verify the 
relevance of the chosen model and the fragment size de- 
pendence of our calculations. 

A least-square fit method is used in order to extract the 
effective parameters from the ab-initio calculations. The 
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conditions imposed for this purpose are that the effective 
model should reproduce 

— the computed excitation energies, 

— the projection of the computed wave-functions within 
the configuration space generated by the magnetic or- 
bitals. 

Let us note that the norm of the later projections give 
us a measure of the model validity. Indeed, in the present 
work, the ab-initio wave- functions (which are expanded 
over several millions of configurations) project over the 
model configuration space based only on the magnetic or- 
bitals (typically of the order of 10 configurations), with a 
norm as large as 0.8. One can thus assume that this model 
space is appropriate to describe the low energy physics of 
the system. 



a) CU2O6 b) Cu 3 8 

ti, Ji, Ae ti, Ji, Ae, V 




t2, J2 



Fig. 1. a) Schematic representation of the computed frag- 
ments, a) two centers b) three centers. The gray circles repre- 
sent oxygen atoms, while the black circles represent the copper 
atoms. 

Let us now address the embedding problem for incom- 
mensurate systems. The usual embedding technique con- 
sists of reproducing the Madelung potential using a set 
of point charges, located at the crystallographic positions 
over a box of at least 15 to 20 A around the fragment |23| . 
These charges are adjusted on the bath borders according 
to an Evjen procedure |24|. In the case of a periodic ionic 
crystal this method insures the nullity of both the system 
charge and dipole moment. These two conditions are nec- 
essary to insure a good representation of the Madelung 
potential. In the present incommensurate case, however, 
the Evjen procedure fails to suppress the dipole moment. 
This is due to the fact that the chain and ladder sub- 
systems are not electrostatically neutral. The relative dis- 
placement of one compared to the other thus induces a 
dipole moment in the chain/ladder direction. In order to 
solve this problem, and to cancel out the dipole moment 
contribution, we adjusted the charges of the outermost 
unit cells of each subsystem, in the c direction, using a 
global scaling factor for each adjusted cell [2HJ. 

The calculations have been done on 11 equivalent frag- 
ments located at 11 successive positions in the chain direc- 
tion. These 11 fragments give a quite good representation 
of the different distortions occurring on the chain sub- 
system. In order to fully represent the whole chain sub- 
system, these 11 sets of results have been extrapolated, 



using Fourier's series analysis, as a function of the crys- 
tallographic fourth coordinate r, which is associated with 
the system incommensurate modulations. 

Let us notice that in a complete crystallographic de- 
scription ^Jjj] , each atom possesses a fourth fractional coor- 
dinate Ti — v\ ■ k = ZiC c /ci, where ri is the atom position 
in the average structure, k = c*c c /q is the modulation 
vector, Zi is the fractional coordinate of the atom in the 
c direction. In the model Hamiltonian used in the present 
work, r corresponds to the fourth coordinate of the chain 
unit-cell copper atom. It is defined except for a constant. 



3 Nature of the holes 

As previously noticed, it has been supposed in the litera- 
ture that the holes (both in the chain and ladder subsys- 
tems) are not located on the copper atoms but rather on 
the surrounding oxygens. This assumption has been done 
by analogy with the high T c copper oxides and has later 
been comforted by NEXAFS experiments as far as 
ladder holes are concerned. We have thus derived from our 
ab-initio calculations the nature and composition of both 
the magnetic (supporting the spins) and hole orbitals. 

The magnetic orbitals have been obtained, from the 
2 centers fragment calculations, as the triplet natural or- 
bitals (eigenfunctions of the one-electron density matrix) 
with an occupation number close to 1. In order to lo- 
cate the hole orbitals we compared two calculations on 
the same fragment with one electron difference. The hole 
orbitals have been extracted from the difference between 
twice the triplet density-matrix and the sum of the two 
doublet ones. The hole orbitals are thus the two eigenfunc- 
tions of the resulting matrix associated with eigenvalues 
close to 1. 

Fig. 2. a) Example of magnetic orbital supporting the spin, 
b) associated hole orbital. 

Figure [3 represents a typical example of both the mag- 
netic (a) and hole (b) orbitals on a site. It clearly appears 
that while the magnetic orbital is essentially supported by 
the 3d ac copper orbital with a derealization tail on the 
surrounding 2p oxygens orbitals, the hole orbital is essen- 
tially supported by the 2p oxygen orbitals with a small 
tail on the 3d ac orbital of the copper atom. The average 
repartition between copper and oxygens — over the 11 
calculations — are as follow : 

— for the magnetic orbital : 67% on the copper and 33% 
on the oxygens, 
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— for the hole orbital : 15% on the copper and 85% on 
the oxygens. 

The fluctuation of the above repartition on the 11 com- 
puted fragments is very small with a standard deviation 
of only 1.3%. This stability of the copper-oxygen com- 
position of the hole and spin supporting orbitals is quite 
noticeable since the electrostatic potential modulation on 
both the copper and oxygen atoms are very strong, with 
an amplitude of 4eV. Finally, we would like to point out 
the stronger oxygen character of the hole orbitals in this 
compound, compared to what was found by Calzado et 
al |27| in the high T c parent compound La^CuO^. Indeed, 
they found only a 50 % oxygen character of the hole sup- 
porting orbital. 

The present results show the existence of different or- 
bitals, respectively supporting the spins and the holes, and 
constitute a direct evidence of the ZR singlets in the sys- 
tem. 



4 Second neighbor t — J + V model 

The parameters of a second neighbor t — J + V model (see 
eq. ^) were extracted from the ab initio calculations on 
the eleven positions of the fragments along the chain. 

Rt-j+v = g(n) th 

i 

a i,cr a i+l,cr + a \+l,a a i,o ) 

i a 

(n) (SiSi+i - 1/4 tiim+i) 

i 

i cr 

~ Wg( T ») ( S i S i+2 - 1/4 niTli+2) 
i 

v( r i) n i n i+i (!) 

i 

where Ti is the fourth coordinate of site i. a\ a , ai, a and rij 
are the usual creation, annihilation and number operators, 
a is the spin quantum number. 

We thus obtained eleven sets of parameters that were 
further fitted as a function of the fourth crystallographic 
coordinate Ti. The fit was done using Fourier series, ac- 
cording to the following expression 

a + a n cos (27rnr - <p n ) (2) 

n 

Only terms with a non negligible contribution to the se- 
ries have been retained. It results that the orbital ener- 
gies, hopping integrals and NN repulsion terms can accu- 
rately be obtained from a unique cosine. Exchange inte- 
grals, however, necessitate two components to be reliably 
reproduced, as can be expected from the quadratic depen- 
dence of the latter to the hopping integrals. 



The results are summarized in table ^ Let us notice 
that only terms with even frequencies have a non negligible 
contribution (n = 2 and n = 4 in eq. |5J). This can be 
interpreted as a doubling of the modulation vector k. The 
fourth coordinate of a unit cell is thus given by 

t = z x 2 c c /ci ~ z x 2 x 7/10 

where z is the fractional coordinate in the c direction of 
the unit cell copper atom. It thus clearly appears that the 
model Hamiltonian presents a 5 unit cells quasi-periodicity. 
This point is in agreement with the neutron scattering ex- 
periments |7lll3| that see a pseudo periodicity of the spin 
arrangement corresponding to five chain unit cells. 
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Table 1. Analytic fit of the t — J + V second neighbor model. 



Figure reports the effective hopping integrals while 
figure 2| reports the effective exchange integrals. Both are 
given as a function of the fourth crystallographic coordi- 
nate t (defined as in the Hamiltonian expression). 
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Fig. 3. NN and NNN effective hopping integrals as a function 
of r. The circles correspond to the two-coppers fragment cal- 
culations, the cross to the three-coppers calculations and the 
lines to the Fourier fits. 

The NN integrals have been evaluated both from the 
two-centers and three-centers fragments. The three-centers 
fragments used in the calculations have been chosen in 
successive positions along the chain, thus, each NN in- 
tegral appears in two successive fragments yielding inde- 
pendent evaluations. One should first notice that the three 
independent evaluations of the the NN hopping, tx, and 
exchange, J\, integrals yield the same values, thus vali- 
dating (i) the second neighbor t — J + V model used in 
the present work and (ii) the strong locality of the interac- 
tions. Indeed, whether other interactions or orbitals would 
have been of importance, whether the integrals screening 
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Fig. 4. NN and NNN effective exchange integrals as a function 
of r. The circles correspond to the two-coppers fragment cal- 
culations, the cross to the three-coppers calculations and the 
lines to the Fourier fits. 



effects would have a larger range than the first coordi- 
nation shell of the magnetic atoms, the two-centers and 
three-centers calculations would have given different eval- 
uations of the integrals so that to compensate the inability 
for the model to reproduce the low-energy local physics. 

As expected, the NNN interactions are of the same 
order of magnitude as the NN ones. In fact, the second 
neighbor hopping integrals are even nearly twice as large 
in amplitude as the first neighbor ones. Indeed, while the 
NN hole orbitals are nearly orthogonal, the NNN hole or- 
bitals strongly overlap (see figure EJa). Another important 
point is the strong modulation of the hopping integrals as 
a function of r. The t-i modulation amplitude is 88 mcV 
and represents a 15% standart deviation around an aver- 
age amplitude of 214meV. 

NN and NNN exchange integrals are, as expected, dif- 
ferent in nature with ferromagnetic NN exchange and an- 
tiferromagnetic NNN one. This is due to the usual well 
known reasons (i) the direct exchange integral, which is 
always ferromagnetic in nature, decreases exponentially 
with the Cu-Cu distance and thus becomes negligible be- 
tween NNN centers, (ii) the super-exchange term, medi- 
ated via the bridging oxygen 2p orbitals and antiferromag- 
netic, is hindered in the NN case due to the destructive 
interaction of the bridging orbitals. It results for the NNN 
exchange a weak, antiferromagnetic integral, in agreement 
with the experimental findings j^j ■ Two points are however 
worth noticing. 

— The experimental evaluation of the intra-dimer exchangi 
Jd is —11 meV, that is somewhat larger than our find- 
ings for NNN exchange. 

— The computed Ji integral strongly vary along the chain 
since its amplitude variation is 55 % of its mean value 



— 7 meV. This strong modulation is however not ob- 
served in neutron scattering experiments 

These two important points are discussed in the following 
sections. 



5 Super-exchange mechanism through a 
double bridge 

The DDCI method used in this work is well known to 
give reliable and very accurate evaluation of the NN ex- 
change integrals. It has been used with a lot of success 
for computing NN exchange integrals in copper |21U27j . 
nickel [2E], vanadium [22J and other oxides. Indeed, it was 
shown that in order to accurately compute magnetic in- 
tegrals, the following effects should be treated [2Hll3()IIHll 

— the multideterminantal and open-shell character of the 
reference wave function, 

— the correlation within the magnetic orbital shell, 

— the screening effects (dynamical polarization) on the 
different configurations of the reference wave function, 

— the effect of the ligand to metal charge transfers con- 
figurations that mediate the interaction, 

— and the screening effects on the latter configurations. 

The characteristic of the DDCI method is to treat all these 
effects in the case of interactions mediated by at most one 
bridging entity. 

In the present case, the NNN interactions are medi- 
ated by two and not one bridging entities, namely the 
two oxygens between the first and second copper and the 
two oxygens between the second and third copper. It re- 
sults a perturbative evaluation of the J2 integral within 
a two bands (magnetic plus bridging entities) extended 
Hubbard model of 

f i f 2 f i f 2 f i f 2 f A f 2 

j = a pd pp _ , L _p£pp_ _ 1fi V> _ 1 fi b pd b pp 

2 A\U d A\A' 2 A\A 2 A' 2 A\A\A< 2 

where t p d is the hopping integral between the copper mag- 
netic orbital and the bridging oxygen 2p orbitals and t pp 
is the hopping between the orbitals of the two bridging 
entities. A\ is the excitation energy between the reference 
configurations and the single ligand-to-metal charge trans- 
fer configurations. A2 and A' 2 are the excitation energies 
of the double charge transfers configurations, Z\ 2 when the 
holes are on different bridging entities and A' 2 when they 
are on the same. Figure [5] pictures the different pertur- 
bative mechanisms yielding the above equation. One sees 
that out of the twenty possible paths, only two of them 
do not involve double ligand-to-metal charge transfer con- 
figurations. Thus, even-though the double charge transfer 
excitations are of higher energy than the single ones (of 
the order of twice as large) , they are much more numerous 
and thus should be responsible for a large part of the J2 ef- 
fective integral, unlike the case of single bridged integrals. 
One should now remember that while the DDCI method 
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Fig. 5. Perturbative mechanism of the NNN effective exchange integrals mediated through two bridging entities. In each config- 
uration, the relative position of the orbitals are pictured in energy scale, that is the two higher orbitals are the magnetic copper 
ones and the two lower orbitals are the 2p bridging entities orbitals. The single ligand-to metal charge transfer configurations 
are frame with dashed lines and the double transfer ts with solid lines. The copper to ligand t p d hopping integrals are pictured 
with solid arrows while the ligand to ligand hoppings, t pp , are with dashed arrows. 



correctly compute the screening effect of the ligand-to- 
metal charge transfer configurations, this is not the case 
for the double charge transfer ones |3D]. It results that the 
excitation energies toward the double charge transfer con- 
figurations (A2 and A' 2 ) are overestimated and the that 
thus the Ji integral is underestimated in the present cal- 
culation. Unfortunately including the screening effects of 
the double charge transfers is technically not possible at 
this time since it would involve triply excited configura- 
tions compared to the reference ones. 

It is important to point out that a similar perturbative 
analysis of the NN integrals as well as the NNN hopping 
integral shows that the contribution of the double charge 
transfers are negligible in these cases. Moreover, noticing 
that the J2 modulations are essentially due do the modu- 
lations of the transfer integrals (t v d and t pp ), the relative 
variations of J2 can be assumed to be correctly reproduced 
in the present calculations. 



6 Second neighbor exchange fluctuations 
along the chains 

Let us now address the question of the NNN exchange 
fluctuations along the chain. As previously mentioned the 
neutron scattering experiments do not see much fluctua- 
tions of the intradimcr effective exchange. In order to com- 
pare our results with the experimental ones one should 
determine which among all the NNN exchange integrals 
along the chain correspond to an intradimcr spin coupling. 

In incommensurate systems, the real space analyses 
present the drawback to never provide a complete descrip- 
tion of the total system. We thus chose to identify the 
dimer localization along the chain not in the real space 
but as a function of the fourth crystallographic coordinate 
t. Even if such an analysis is not as intuitive, it presents 
the advantage to provide a complete description of the 
system. 

We will thus analyze the system in order to derive a 
partition into physically pertinents blocks and to deter- 
mine (i) the number of different type of blocks, (ii) their 



relative arrangement and (iii) their rate of existence in 
the incommensurate structure. For this purpose it is use- 
ful to start with a chosen type of block, expected to be 
largely represented in the system. In the present case we 
chose the dimeric units. Then, the model Hamiltonian is 
used to determine the apparition of this type of blocks as 
a function of the fourth component. The analysis of the 
remaining sites allows to determine the other types of sig- 
nificant blocks present in the system. Finally, the relative 
arrangement of the types of blocks can be studied. 

The dimeric units are composed of five sites, the three 
sites of the dimer plus the two ZR singlets of the neighbor- 
ing sites (see figure EJ. Let us use the fourth coordinate 
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dimer blocks (85 %) 



free spin blocks (15 %) 



12 3 



Fig. 6. The different types of spin blocks appearing in the 
chains. Values in parentheses refer to the percentage of the 
chain corresponding to the different types of blocks. 



of the first site (rx) as the block reference. If one plots 
the energy of five consecutive sites as a function of t\ , a 
dimeric unit will thus be obtained when the energy of the 
second and fourth sites are below the Fermi level, while 
the energy of the first, third and fifth sites are above the 
Fermi level. Figure shows the orbital energy curves of 
five consecutive sites as a function of T\. Let us notice 
that, since five sites are represented by a unique value of 
t, only one fifth of the [0, 1] range is necessary to represent 
the whole chain. The values of t\ for which the five con- 
secutive sites form a dimer are represented in gray. They 
span two ranges of width 0.085. The proportion of the 
chain occupied by dimers thus corresponds to five time 
the above ranges, that is 85% of the whole chain. The 
number of dimers can then be evaluated to 1.70 per f.u., 
to be compared to 1.47 dimer per f.u. deduced from mag- 
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Fig. 7. Orbital energy of five consecutive sites as a function 
of the fourth coordinate n of the first one. The horizontal line 
represent the Fermi level. 
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Fig. 8. Orbital energy of block sites for the two different types 
of block given as a function of the fourth coordinate n of the 
first site of each block. The horizontal line represent the Fermi 
level. 



netic susceptibility measurements Let us now analyze 
the composition of the remaining 15% of the chain. For 
this purpose we will determine the number of sites sep- 
arating two consecutive dimers. If t\ is referencing the 
first dimer (that is n € [0.068; 0.153] U [0.568; 0.653]) 
then the reference of the second dimer is given by r{ = 
ri+n(ri)c c /c;, where n{j\) is the smaller integer such that 
t{ e [0.068; 0.153] U [0.568; 0.653]. At this stage three 
cases can occur. 

1. case, ti(ti) < 5. In this case, the two successive dimers 
overlap, in other words there exist physical entities 
larger than the dimers. 

2. case, ti(ti) = 5. In this case, the two dimers are strictly 
consecutive along the chain. 

3. case. n{ri) > 5. In this case, the two successive dimers 
are separated by one or several other type of blocks of 
total length : n(n) — 5. 

In the present system, only n = 5 and n = 8 occurs. As 
stated above the former corresponds to consecutive dimers 
and the latter to dimers separated by three sites blocks. 
As for the dimers, we will reference these blocks by the 
fourth coordinate of their first site, that is t\ + 5 c c /q . The 
three sites blocks thus span two ranges : [0.042 ; 0.068] and 
[0.542; 0.568]. Figure |H1 represents the three sites blocks 
and the dimer ranges for n < 0.5. In each region, the 
orbital energy of each site of the block has been repre- 
sented. One sees immediately that the three sites blocks 
are formed by one spin surrounded by two ZR singlets. 
Such a configuration can be associated with free spins 
since the nearest neighbor spin is two ZR singlets afar. 
Let us notice that three times the free spins ranges yields 
0.15, that is the total missing part of the chain. The num- 
ber of free spins can be now easily evaluated to 0.5 per 
f.u., to be compared to 0.55 free spins per f.u. obtained 
from the magnetic susceptibility experiments. 

A further analysis of the figure[7|shows that the dimers 
are arranged in clusters of three or four dimers separated 
by a free spin. It can be evaluated that 54% of the dimers 
form three-dimers clusters while 46% form four-dimers 
clusters. 



Let us now go back to the intradimer integrals. The 
regions of r corresponding to dimers are reported in gray 
on figure El The intradimer J2 regions correspond to J2 




Fig. 9. NNN effective exchange integrals as a function of r. 
The regions in gray correspond to intradimer effective exchange 
integrals. Let us note that the r values in the present figure is 
associated with the T2 value of the second site of the 5 centers 
dimeric units. This results in a shift in the r values compared 
to figure [7| 



values between —8.1 and — 8.8meV. 

Our calculations thus yield values of the intradimer ex- 
change that are quite homogeneous since their fluctuation 
range is smaller than 0.7 meV, that is less than 8% of their 
nominal values, in good agreement with the neutron scat- 
tering experiment (the variation range is within the ex- 
perimental error bars) . The other interesting point is that 
the Ji values corresponding to intradimer exchanges are 
the largest ones in absolute values. The underestimation 
due to our method is thus quite limited since our average 
value is —8.5 meV to be compared to the experimental 
values ranging from 11 and 12 meV. 



7 Conclusion 

To summarize the present results, we have determined a 
second neighbor t — J + V model for the incommensu- 
rate chain subsystem of the Stx^Cu-iaQ^x compound. The 
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model parameters have been determined using accurate ab 
initio calculations on a series of embedded clusters along 
the chain. In order to obtain a complete model as a func- 
tion of the incommensurate modulations, the ab initio re- 
sults have been extrapolated using a Fourier analysis. The 
resulting model is thus independent of any periodic ap- 
proximation and is given as a continuous function of the 
fourth crystallographic coordinate r, which describes the 
incommensurate modulations along the chain. 

It is noticeable that, unlike what is currently assumed 
in the literature, the various parameters of the model, ex- 
cept for the first neighbor bi-electronic repulsion, strongly 
vary as a function of the structural incommensurate mod- 
ulations. In fact, these variations are so large that they 
determine the physics of the system. In particular, the 
orbital energies vary over a surprisingly large range, and 
thus dominate the low energy physics through a strong lo- 
calization of the electron (resp. holes) over the low (resp. 
high) energy sites. 

The analysis of the model as a function of the fourth 
crystallographic coordinate r allowed us to show that at 
low temperature the chain ground state can be entirely 
described only by second-neighbor dimers and free spins 
if no hole transfer to the ladder is assumed. It is worth 
to point out that our estimation of the free spin number 
per formula unit, 0.5, is under these filling hypothesis in 
very close agreement with the experimental value of 0.55. 
Further analysis of our results showed that the dimers are 
arranged in clusters of three or four units separated by a 
free spin. The intradimer exchanges are found to be an- 
tiferromagnetic, weakly dependent of the dimer position 
along the chain and in reasonable agreement with the ex- 
perimental values, even if a little underestimated. It is 
amazing to realize that while the strong modulation of 
the model parameters are responsible of the low energy 
properties of this system and in particular of the elec- 
tron / hole localization and the resulting formation of the 
dimers, the parameters range actually seen in magnetic 
experiments corresponds to a nearly homogeneous dimer 
system. 

Wave functions analysis exhibited holes localization 
not on the copper atoms but rather on the surrounding 
oxygen 2p orbitals, thus confirming the hypothesis of the 
presence of Zhang-Rice singlets. Three types of Zhang- 
Rice singlets can be identified in our calculations, namely 
the intra-dimer ones, and two types of inter-dimer ones : 
those neighboring dimers and those neighboring free spins. 
This result is to be put in perspective with the copper 
NMR experiments [4] that sees a splitting of the inter- 
dimer Zhang-Rice singlets signal at low temperatures. It 
would be interesting to quantify the relative weight of the 
two signals in order to check our predicted ratio of 3.4. 

Acknowledgment : the authors thank Dr. J. Etril- 
lard for providing us with the neutron crystallographic 
structures as well as for helpful discussions, Dr. D. Maynau 
for providing us with the CASDI suite of programs, Dr. P. 
Sciau for introducing us with the four-dimensional crys- 
tallographic conventions. The present calculations were 



done using the IDRIS computational center facilities un- 
der project number 1104. 



References 

1. M. McCarron et al., Mater Res. Bull. 23, 1355 (1988) 

2. S.A. Carter et al, Phys. Rev. Lett. 77, 1378 (1996). 

3. R.S. Eccleston, M. Azuma and M. Takano, Phys. Rev. 
B 53, 14721 (1996) ; M. Matsuda et al, Phys. Rev. B 
54, 12199 (1996). 

4. K. Magishi et al, Phys. Rev. B 57, 11533 (1998). 

5. S. Tsuji et al, J. Phys. Soc. Jap. 65, 3474 (1996) ; 
K.Kumagai et al, Phys. Rev. Lett. 78, 1992 (1997). 

6. M. Takigawa et al, Phys. Rev. B 57, 1124 (1998). 

7. R.S. Eccleston et al, Phys. Rev. Lett. 81, 1702 (1998). 

8. U. Ammerahl et al, Phys. Rev. B 62, 8630 (2000). 

9. M. Kato, K. Shiota and Y Koike, Physica C 258, 284 

(1996) . 

10. F. C. Zhang and T. M. Rice, Phys. Rev. B 37, 3759 
(1988). 

11. N. Niicker et al, Phys. Rev. B 62, 14384 (2000). 

12. D.E. Cox, I. Iglesiais, K. Hirota, G. Shirane, M. Mat- 
suda, N. Motoyama, H. Eisaki and S. Uchida, Phys. Rev. 
B 57, 10750 (1998). 

13. LP. Regnault et al, Phys. Rev. B 59, 1055 (1999) ; M. 
Matsuda et al, Phys. Rev. B 59, 1060 (1999). 

14. C. Hess, H. El Haes, B. Biichner, U. Ammerahl, M. 
Hiicker a,d A. Revcholevschi, Phys. Rev. Letters 93, 
27005 (2004). 

15. C. K. Majumdar and D. K. Gosh, J. Math. Phys 10 
1399 (1969). 

16. F. D. M. Haldane, Phys. Rev. B 25 4925 (1982). 

17. A. Gelle and M.-B. Lepetit, Phys. Rev. Letters 92, 
236402 (2004). 

18. M.-B. Lepetit, Recent Research Develepments in Quan- 
tum Chemistry, p. 143, Vol. 3, Transworld research Net- 
work (2002) (or |cond-m at/0502268ll. 

19. N. W. Winter, R. M. Pitzer and D. K. Temple, J. Chem. 
Phys. 86, 3549 (1987). 

20. J. Miralles, J. P. Daudey and R. Caballol, Chem. Phys. 
Lett. 198, 555 (1992) ; V. M. Garcia et al, Chem. Phys. 
Lett. 238, 222 (1995) ; V. M. Garcia, M. Reguero and 
R. Caballol, Theor. Chem. Acc. 98, 50 (1997). 

21. D. Munoz, F. Illas and I. de PR. Moreira, Phys. Rev. 
Letters 84, 1579 (2000) ; D. Munoz, I. de PR. Moreira 
and F. Illas, Phys. Rev.B 65, 224521 (2002). 

22. N. Suaud and M.-B. Lepetit, Phys. Rev. B 62 402 
(2000) ; N. Suaud and M.-B. Lepetit, Phys. Rev. Letters 
88, 056405 (2002). 

23. M. Isobe et al, Phys. Rev. B 62, 11667 (2000) ; J. Etril- 
lard et al, Physica C, in press. 

24. H. M. Evjen, Phys. Rev. 39, 675 (1932). 

25. A. Gelle and M.-B. Lepetit, unpublished. 

26. P.M. de Wolff, Acta Cryst. A 30, 777 (1974). 

27. C. J. Calzadoa, J. F. Sanz and J. P. Malrieu, J. Chem. 
Phys. 112, 5158 (2000). 

28. I. De P. R. Moreira and F. Illas, Phys. Rev. B 55, 4129 

(1997) ; C. de Graaf, F. Illas, R. Broer, and W.C. Nieuw- 
poort, J. Chem. Phys. 106, 3287 (1997). 

29. P. de Loth, P. Cassoux, J. P. Daudey et J. P. Malrieu, 
J. Am. Chem. Soc, 103, 4007 (1981). 



Alain Gelle, Marie-Bernadette Lepetit: SK4CM24O41 : a complete model for the chain sub-system 



9 



30. A. Gelle, M. L. Munzarova, MB Lepetit and F. Illas, 
Phys. Rev. B 68, 125103 (2003). 

31. C. J. Calzado, J. F. Sanz and J. P. Malrieu, J. Chem. 
Phys. 112, 5158 (2002). 

32. MB Lepetit, Recent Research Developments in Quan- 
tum Chemistry 3, p. 143, Transword Research Network 
(2002). 



